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Abstract. We consider an integrable three-dimensional system of ordinary differential 
equations introduced by S.V. Kovalevskaya in a letter to G. Mittag-Leffler. We prove 
its isomorphism with the three-dimensional Euler top, and propose two integrable dis- 
cretizations for it. Then we present an integrable generalization of the Kovalevskaya 
system, and study the problem of integrable discretization for this generalized system. 



Or 



1. Introduction 

The problem of finding integrable discretizations for integrable systems of ordinary 
differential equations was studied a lot in recent years, see [llj. Depending on the in- 
tegrability attributes one wants to respect, different approaches can be proposed and 
sometimes different results can be achieved for a given integrable system. One common 
feature of almost all integrable discretizations known to date is that equations of motion 
of these maps appear to be rather nontrivial deformations of the original differential equa- 
tions. These deformations can only be found after running a sometimes involved machine 
of some systematic approach, or, alternatively, by a guesswork based on an extensive 
experience and expert knowledge. 

We are aware of only one approach for which the situation is different. This approach, 
due to R. Hirota and K. Kimura [3], is applicable to any vector field with a quadratic 
dependence on coordinates. It consists in replacing the derivatives by the first differences 
and the quadratic terms in the vector fields by the corresponding bilinear terms (with re- 
spect to the values of coordinates in two subsequent instances of the discrete time). Thus, 
a derivation of a Hirota-Kimura discretization for a system is quite straightforward and 
algorithmic. However, there are no theoretical reasons for Hirota-Kimura discretizations 
to preserve integrability. Nevertheless, very remarkably, it turns out that very often they 
do! An extensive experimental and theoretical study of this phenomenon is contained 
in our recent works [UEHE]- These papers contain a big collection of novel integrable 
birational maps, whose integrability is expected to be related to some deep phenomena 
in algebraic geometry. 

The goal of the present work is to add several new instances to this collection, with the 
hope that this will facilitate the development of the still lacking theoretical explanation 
of these exciting observations. 



E-mail: petrera@math.tu-berlin.de, suris@math.tu-berlin.de. 
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2. KOVALEVSKAYA SYSTEM AND ITS RELATION TO THE EULER TOP 

In a letter to G. Mittag-Leffler from December 1884 [5_> p. 80-82], S. Kovalevskaya 
discussed a class of systems of ordinary differential equations with quadratic vector fields, 

3 

V% = yi^2 a ijVj, i = 1,2,3, (1) 
i=i 

where are parameters. She mentioned that under the condition 012023031 = 013032021, 
a general solution of such a system depends on three arbitrary parameters, and its only 
singularities in the finite part of the complex plane of the independent variable are first 
order poles. In the modern theory, this property is called "Painleve property" and is 
associated (somewhat vaguely) with integrability. As a particular case of this class, she 
mentioned the following system: 

yi = 2/i(— 2/1 + 2/2 + 2/3), 

yi = 2/2(— 2/2 + 1/3 + yi), (2) 
h = 2/3(— 2/3 + 2/1 + 2/2)- 

This system admits three conserved quantities: 

^23 = 2/l(2/2 — 2/3), ^31 = 2/2 (2/3 — 2/1 ) , Kl2 = 2/3(2/1 -2/2)- (3) 

They are linearly dependent, K23 + -K31 + K\ 2 = 0, and any two of them are functionally 
independent. In her letter, S. Kovalevskaya stated that system (T5]) can be integrated in 
terms of elliptic functions and suggested that the study of the class of systems (pP) possess- 
ing the Painleve property might lead to a better understanding of general systems with 
quadratic vector fields. Although this opinion is probably too optimistic, it is definitely 
useful to have a closer look at system ([2]) and its integrable generalizations. 

What S. Kovalevskaya did not mention is a simple (and even two- fold) relation of 
system fl2j) to a much more famous integrable system, the Euler top (see, e.g., [TU]): 

xi = x 2 x 3 , 

x 2 = X3X1, (4) 

£3 = ^X 2 . 

One can bring the latter system to the form Xi = otiXjXk, involving arbitrary real parame- 
ters Oj, by a simple scaling transformation x« t— > PiXi, where, however, some of may be 
imaginary, depending on the signs of ctj. Here and below we often use k) to denote 
a cyclic permutation of (1,2,3). 

System (J5J) admits three integrals of motion 

Eij X ^ Xj , 

only two of which are functionally independent due to E 2 3 + E31 + E 12 = 0. It is easy to 
verify that both the following changes of variables, a linear one, 

Vi = gfe + x k) & x i = ~yi + Vj + 2/fc, (5) 

and a nonlinear one, 

Vi = & Xi = VWE, (6) 
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map (jl]) to (j2]). Indeed, for the change of variables (j^J) we compute: 

1 1 

yi = g w + = 2 X ^ X ^' + Xk ^ = y^y* + yj + ^ fc )> 

while for the change of variables ([6]) the corresponding computation gives: 

1 / . <y . hp , iy . iy . y , iy . y . >y . y , 

— = — + = — + — ~ - J — = -Vi + Vi + y*- 

i/i Xj Xk Xi Xj Xk Xi 

One easily sees that the conserved quantities Ky of (T5]) are obtained from the correspond- 
ing conserved quantities Eij of (jl]) by either of the changes of variables. 



3. TWO INTEGRABLE DISCRETIZATIONS OF THE KOVALEVSKAYA SYSTEM 

Any system of ordinary differential equations with a quadratic vector field can be 
discretized via the Hirota-Kimura (HK) approach, which consists of the replacement of 
the derivatives Xj by the differences (xj — Xj)/e and the quadratic expressions like XjXk 
or x 2 by their respective bilinear counterparts, XjX k + XjXk, resp. %XjXj. For a function 
x : eL — > R, we use the abbreviation x for x(t) and x for x(t + e), where a is the discrete 
time-step. We refer to jUEHS] for our recent investigations on HK-type discretizations. 

The HK-discretization of the Euler top (jlj) is described by the following discrete equa- 
tions of motion [3]: 

x\ — x\ = e(x 2 x 3 + x 2 x 3 ), 

x 2 - x 2 = e(x 3 x 1 + x 3 xi), (7) 
x 3 - x 3 = e(x 1 x 2 + x\x 2 ). 

These equations define a birational map (xi,X2,Xs) = f(xi,x 2 ,x 3 ,e), reversible because 
of e) = f(x, —e). Its explicit form is: 

*£ 2 I 2 £ Xi j Xi I £ Xi ^ ^ Xj ^ I Xj A I xi.) 



s 2 (x\ + x\ + x\) — Is^xxx-iXz 



As proven in [8], this map is bi-Hamiltonian and completely integrable. Its orbits are 
given in terms of elliptic functions, and it admits conserved quantities 



^v>2 ,y 2 

W (r\ — —" 



ran \~ i -, 09) 
E Z Xj 



with arbitrary m,n,j = 1, 2, 3. (This notation is incomplete as it does not refer to j; only 
the indices surving in the continuous limit e — > are explicitly referred to.) Among the 
functions E mn (e), there are two functionally independent ones. Moreover, the map (JED 
has an invariant volume form 

dxi A dx 2 A dx 3 



OJ 



<j>{x) 



where <fr(x) is any of the functions (1 — e 2 x|) 2 . Note that the ratio of any two such 
functions is a conserved quantity. 

It turns out that the HK-discrete Euler top ([7]) can be transformed into birational 
integrable discretizations of the Kovalevskaya system (J2J) by any of the changes of variables 
(EJ) and However, the resulting maps are different. 
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Linear changes of variables preserve the class of HK-discretizations. Therefore, the 
change of variables (0) maps the HK-discrete Euler top onto the HK-discrete Kovalevskaya 
system: 

yi-yi = £{-2ym + ym + ym + ym + ym), 
yi-yi = e{-2y 2 y 2 + mys + ym + fayi + 2/22/1), (9) 

m - 2/3 = e(-2y 3 y 3 + y 3 y x + y 3 y x + y 3 y 2 + y 3 y 2 ). 
The properties of the HK-discrete Euler top are then translated as follows. 

Theorem 1. 

(1) The map (TJ|) admits the following conserved quantities: 

K mn (E) = l _ e 2( j/ ^Ty j , + ytf • 

where K mn , with arbitrary and distinct m,n — 1, 2, 3 are defined in ^ and k) 
is any permutation of (1,2,3). Among functions [TU) only two are functionally 
independent. 

(2) The map admits an invariant volume form 

_ dyi A dy 2 A dy 3 

where if)(y) is any of the functions (1 — e 2 (%)i — yj + yk) 2 ) 2 ■ 

Now we apply the nonlinear change of variables ([6]) to the HK-discretization of the 
Euler top (J7J). A remarkable feature of the map (JHJ), written as 



1 + 2e(xjX k /xi) + s\-xj + xj + x{) 



Xi 1 — e 2 ix\ + x\ + x\) — 2e 3 xix 2 x 3 

is that the r.h.s of the latter formula depends rationally on combinations of Xi of the form 
XjXk/xi = yi, x 2 = yjyk, and x\x 2 x 3 = yiy 2 y 3 , which do not involve square roots Jy]. As 
a consequence, the pull-back of the HK-discrete Euler top under change of variables (jBJ) 
is birational: 

~ _ y ( 1 + 2e yi + g2 (^i + ViVk ~ VkVi)) (l + 2ey k + e 2 (yjy k + y k yi - ViVj)) ^ 
1 (l + 2ey { , + e 2 (y k yi + yiyj - y^k)) (l - £ 2 {ym + y 2 y 3 + y 3 y x ) - 26^^^) 



Before translating properties of the HK-discrete Euler top into properties of map (ITT1) . 
we mention a result of our recent work [9], according to which the HK-discretization of 
the Euler top is the second iterate of another map which is nothing but the cosine law 
for spherical triangles in § 2 . Indeed, the cosine law can be written as 

x i ~t~ x j x k 

Xi 



where Xi = cosaj, Xi = cos£j, with a, and li being angles and sides, respectively, of a 
spherical triangle. The small parameter e can be inserted by a scaling transformation 
which makes out of the latter map a "new" discretization of the Euler 

top: 

— . _ X i + ex j x k /-,„v 
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The second iterate of (jlzj) . where defined, is exactly the birational map (J8J). 
Now, performing the nonlinear change of variables ([6]) in ffT2"l) . we easily get: 



(l + eyj .)(l + eyfc ) 
V * ~ V \l + e Vi )(l - e* yj y k y 



which, remarkably, is birational. The second iterate of ffTBj) is, of course, the map (TTTT) . 
One can say that the latter map admits a birational "square root" map. 

The properties of the HK-discrete Euler top are translated under (jH]) as follows. 

Theorem 2. 

(1) The map / fig)) and its second iterate ITU]) admit the following conserved quantities: 

K mn {e) = = K T , (14) 

where (m,n) and are two arbitrary pairs of distinct numbers from {1,2,3}, 
and K mn are defined in (E|). Among functions ( [1^] ) tnere are two functionally in- 
dependent ones. 

(2) TTie map / fig)) ana 1 zfo second iterate / TIT]) admit an invariant volume form 

_ dyi A cfo/ 2 A dy 3 

where ip(y) is any of the functions (1 — e 2 yiyj) 2 . 

One sees immediately that the ratio of any two functions K mn from flS]) serves as an 
integral of motion for the both discretizations (JHJ) and (Tl3|) . Moreover, they also share an 
invariant volume form with an e- independent density 



4. Generalized Kovalevskaya system 

The main subject of investigation in the present paper is the following system with a 
quadratic vector field, which is a generalization of the Kovalevskaya system (j2J): 

yi = yi{-2 yi + s), l<i<N, (15) 

with 

N 
3=1 

System (TT5|) was proposed in pQ p. 55-56], where, however, only one new result was 
established for this system, namely existence of three polynomial conserved quantities for 
N = A: 

Pi = (yx-y2)(y3-yi), P2 = (vi-y3)(v2-Vi), Ps = (yi- yi) fa-ys)- (17) 

Actually, only two of them are functionally independent, because of P\ — P2 + -P3 = 0. 
One can relate system ([TBI) to an elegant generalization of the Euler top proposed in [2] : 

ii = Y[xj, l<i<N. (18) 
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System ()18|) preserves the phase volume (since its vector field is obviously divergence-free) 
and admits a set of conserved quantities 

Ea = x] - 4 (19) 

Since one can find N — 1 functionally independent conserved quantities among Eij, system 
( ITS]) is superintegrable. Moreover, it can be explicitly integrated as follows: use the 
integrals to express N — 1 variables, say Xj with 2 < j < N, in terms of one variable, say 
xi, leading to a differential equation for x\\ 

N 

(io a =n(^+%)- 

i=2 

Thus, integration is performed through the inversion of a hyperelliptic integral. 

Systems (|15[) and (ITS]) are related by a change of variables which is a direct general- 
ization of the nonlinear change of variables 



vi = ^U x i * ^=£(Ek) ■ ( 2 °) 

* 4 \i=i / 

This is shown by a direct computation: 

- = -- + J2 — = -y* + = + s - 

The properties of system (fl8|) are translated into properties of system (TToT) as follows. 

Theorem 3. 

(1) The functions 

a*) <2i) 

are conserved quantities of system / T73]) . Among functions (GOP, t/iere are iV — 1 
functionally independent ones. 

(2) System < T73]) admits an invariant volume form 

N-3 

dy x A • • • A dyjv , , /-^r \ iV " 2 . . 

w= — , 0U/)=M]_2/fc • (22) 



vfe=l 



0(y) 

Proof. The first statement is an immediate corollary of ( 119]) and ( 1201) . The second follows 
from an easily verified formula 

**(s) = (n*j in- 

Thus, system (fl5]) is superintegrable for any iV. In particular, for N = 4 one has three 
independent integrals among 



V 
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where (i,j,k,£) is any permutation of (1,2,3,4). Integrals fTT7|) are obtained from those 
as follows: Pi = -^12-^34, Pi = K 13 K 2 a, P3 = K 14 K 23 - 

We remark that system ([151) admits a further integrable generalization: 



i)i = yii-ayi + s) 



Ki< N, 



(24) 



with an arbitrary real a 7^ N. A direct computation shows that (1241 has N—l functionally 
independent conserved quantities, contained in the set of functions 



Ki. 




(25) 



and an invariant volume form 

dyi A • 



A dy N 



N 



N+l-2a 
N-a 



<P{y) = [Ylvk 



vfc=l 



Moreover, one can establish the existence of N — 2 independent integrals of motion for a 
still more general system, namely for (1241 with any symmetric function s, not necessarily 
given by ([TBI) . Indeed, setting = (y, — y^/fyiyj), we immediately derive from 



= -sH, 



Vi Vj \Vi Vi. 

Therefore, the functions Hij/Hkt = Kij/K^ are conserved quantities irrespectively of s. 

5. HlROTA-KlMURA DISCRETIZATION OF THE GENERALIZED KOVALEVSKAYA SYSTEM 



Unlike in the case N = 3, the relation to the top-like system (Tl8|) does not help to 
find integrable discretizations for (Tl5|) . since we are not aware of any decent integrable 
discretization of the former system. On the contrary, since system (Tl5|) has a quadratic 
vector field, we are in a position to apply a general Hirota-Kimura discretization scheme 
to it. The resulting discrete equations of motion read: 



Vi -y% = e (-^ViVi + yiS + yis) 



1< i < N, 



(26) 



where s is defined as in (fT6|) . As usual for HK-discretizations, equations of motion (T26l) 
can be solved for y, yielding the rational map 

V = f(y, e) = A-\y, e)y, V = (yi, ■ ■ ■ , Z/tv) T , 



(27) 



with 



A{y,e) 



-zy2 



-eyi 

1 -e(-3j/ 2 + s) 



-eyi 



-ey N -i -£Vn-i ■ ■ ■ -£Vn-i 

\ ~ey N -eyn ■■■ 1 - e(-3y N + s) J 

This map is actually birational, due to the reversibility property / _1 (?/,£:) = f(y, 

Theorem 4. Map (2l\ ) can be represented in the following explicit form: 

1 Vi 



S(y,e) d{' 



1< % < N, 



(28) 
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where 

di = l-e(-4yi + s), (29) 

and 



N 



S(y,e) = l-eJ2¥- (30) 



j= i J 



Proof. The matrix A(y,e) may be written as 



N, 



A(y } e) = D — eye , D = diag (d\, . . . , d 
where e = (1, . . . , 1) T . Therefore, 

y = A- 1 (y J e)y = (l - eD~ y T ) D^y 

l-e(D x y,e)) 

(D- l y,e) \ D- l y 



V l-e(D x y,e)) 



l-e(D-iy,e)J S(y } e) 
where 5 = 1— e(D~ l y, e). This coincides with f l28|) . 

Theorem 5. The following formulas hold true: 



Proof. From (1281) and (130)1 we derive: 

iV 

= 1 - e^2S(y,e)yj, 

3=1 

which gives S(y,e)(l + es) = 1. Similarly for the second formula. 

Theorem 6. 

(1) The functions 

Ky _y%- yj ykyi 



l<iJ,k,£<N, (32) 
K M yiyj yk - ye 

are conserved quantities of the map l[2h)) . Among them, there are N — 2 function- 
ally independent ones. 

(2) The map l[2h)) admits an invariant volume form 

dyi A ■ ■ • A dyx 

where ip(y) is any of the functions 



N-l / N ^ 2 



= III ] < 33 > 



\k=l 
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Proof. (1) Equations of motion fl26l) are equivalent to 

1 + ss 1 — es 



(34) 



As. 

Vi Vi 

Therefore, 

(1 + eS) (I - I) = (1 - ea ) (i-I 

or 

Vi ~ Vj ViVj 1-es 
ViVj Vi ~Vj 1 + es' 
The r.h.s. of this formula does not depend on indices which proves the first statement 
of the Theorem. 

(2) We start with the following relation, which holds true for any HK-discretization: 

Determinants on the r.h.s are easily computed: 

det A(y,e) = det(D) det(l - eD~ 1 ye T ) 

N j iV 

fe=i ; fe=i ^ 

In this computation we used the fact that det(l — eD~ 1 ye T ) — 1 — e {D~ x y, e) = e), 
as well as fl28|) . Upon using the first of equations in fl3~Tl) . we find: 

JV 

detA(y,e) = (l + e^ N - 1 T\ 1 ^. 

Vk 
k=l yK 

Similarly, we have: 

N ~ 

det^(y,- e ) = (l- ea )^- i n-- 

Li y« 

As a consequence, we find: 



With the help of ( I34p this can be represented as 



iV-l / iV 



detff \ = (*^JW \ " TT^ 

v<%/ V Vi-yjJ \t=i yk , 

which is equivalent to the second statement. ■ 

Remark 1. Note that the ratio of any two of distinct functions ip(y) in (133]) is an integral 
of motion, according to statement (1) of Theorem [6j Note also that functions ip{y) do not 
depend on e. Therefore, they must serve as invariant volume densities for the continuous 



■ij \ 

invariant volume density for the flow ([TBll . given respectively in (T2T]) and (122]) . □ 



flow (fl5]) . as well. Indeed, ip(y) = K { - (y)4>(y), where and are integrals and an 
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dfc=l 



Remark 2. Statement (1) of Theorem [6] holds true, exactly in the same form and with 
the same proof, for the HK-discretization of the system ([24"j) with arbitrary a. Moreover, 
the specific form of the function s is also irrelevant for the claim and for the proof, as 
long as it is a symmetric function of the coordinates. □ 

Now the question is whether the map (1261) admits the last, (N — l)-th, integral of 
motion. Conserved quantities (1211) from Theorem [3] suggest to look for it in the form 

K i:j (e) = K i:i $(y,e), 

with some <&(y,e) which is an even function of e. It follows from (j34j) that the necessary 
and sufficient condition for this is given by the functional equation 

${y,e) l-es 

It turns out that a solution of the latter equation can be always found as a formal series 

oo 

= i + 5> 2 W?/). 

k=l 

However, it seems that this series converges (and defines a genuine function) only in the 
cases N = 3,4. Thus, only in these cases can we claim that the discrete system (j2"6"|) admits 
the same number of conserved quantities (2 and 3, respectively) as the corresponding 
continuous system. 

In the case N = 3, the existence of two independent integrals of motion has been 
already given in Theorem [1] as a consequence of the integrability of the HK-discretization 
of the Euler system (151) . We found that the map (126]) with N = 3 admits the following 
conserved quantities: 

K mn (e) = K mn $(y,e), *(y,s) = { _ ^- 1- + . 

Here K mn are defined in (j3J), while (i,j,k) is any permutation of (1,2,3). This set of 
integrals contains two functionally independent ones. 
In the case N = 4, we have the following result. 

Theorem 7. The map (TJj2J) with N = 4 admits the following conserved quantities: 

K mn (e) = K mn e), $(y, e) = 1 (35) 

V 1 - e%< + yj -y k - yt) 

Here K mn , with arbitrary m,n — 1,2,3,4, are defined in ( [Hp , while (i,j,k,£) is any 
permutation of (1,2, 3, 4). Among functions there are three functionally independent 
ones. 



Proof. We have to prove that 

'Vm-Vn VrnVn \ 2 Vi^VzVa _ 1 - £ 2 (yi + Vj ~ Vk ~ Vi) 2 



VmVn Vm-Vn) VmVzVA 1 " £ 2 {Vi + Vj - Vk ~ Vl) 2 

As a direct consequence of (J2"8j) and fTSU]) . we have: 



(36) 



~ ~ 1 — £s Vi ~ Vj 

Vi - Vj = ^7 V a a ■ 37 
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This formula, together with ( 128)1 . yields the following expression for the l.h.s. of (J36J): 

Vra-Vn VmVn \ ^ Vl&VzVA (1 ~ Es) 2 I 



VrnVn Vm ~ Vn J 2/12/22/32/4 S 2 did 2 d 3 d 4 ' 

To find a proper expression for the r.h.s. of (1361) . we start with the following formula for 

^ 4 \di + d 2 + d 3 + d±) 

It is obtained upon taking into account ( 130)1 and expressing y$ = (d$ — (1 — es))/(4e) by 
virtue of ( 125)1 . Now, using ( 157)1 to express the differences yi — yk and yj—yi and expressions 
y% — yk = (di — djfc)/(4e) following from (125)1 . we have: 



'1 - es) 2 (\ 1 1 1 N 2 



16S 12 \dj dfc dj da 

Upon using (155)1 . we find: 

Using again expressions y-i — yu = {di — dfc)/(4e) and relation di + d 2 + d 3 + d 4 = 4, we 
get: 

1 - e 2 (yi + % - y k - = 1 - — (di + dj -d k - d e ) 2 = ~(dj + d j )(d fc + d e ), 
Thus, we find: 

1 - e 2 (y i + y j -y k -y E ) 2 _ (1 - es) 2 1 



1 - e 2 {yi + yj -y k - ye) 2 S 2 did 2 d 3 d 4 : 

which coincides with 



6. An alternative DISCRETIZATION of THE GENERALIZED KOVALEVSKAYA SYSTEM 

We now present a discretization of the generalized Kovalevskaya system flT5)) which 
generalizes ( 1T5)) . It is different from the HK-discretization (|26p . but shares with it N — 2 
integrals of motion and the invariant measure form: 

Vi Hi 1 



l-eyi 1 + eyi R(y, e) 

where 



1 < i < N, (40) 



N 



R( y ,e) = l-eJ2j 



j= i ' ' ^ 



This defines a rational map y = f(y,s). Notice that 

-i 



/ N \ 1 , d\ N 

R(y,e) = n( i+ ^) ( i -^)ri( i +^) 

\j=l / \ / jr = l 



\.y= 

v \ _1 

11(1 + 6^)1 D(y,e), 
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where 

N 

D(y,e) = l-Y,z k (k-l)e k (y), 

k=2 

and e k (y) is the elementary symmetric polynomial of degree k in N variables y: 

i<ii<-<ik<N 

Equations of motion ( 140]) immediately yield also 

Vi = * v l<i<N, (41) 

l + eyi Ri{y,e) 



where 



Ri(y,e) = 



y.i 



+ £%•)) A(2/,e), 



with 



A(2/,e) = D{y,e)\y t =o = 1 -^£ fc (/c- l)e fc (yi, . . . . . . ,y N ), 

k=2 

where y{ means that the variable yi is omitted from the list of arguments. 
For N = 3 map (j4"Ij) looks as follows: 

~ = (1 + gyj)(l + ey k ) 
l {^ + ey^{l-e 2 yjykV 

where k) is any permutation of (1,2,3). This is exactly the map (|T3|) obtained from 
the spherical cosine law. 
For N = 4 we get: 

~ = {l + £yj){l + £y k )(l + £y e ) 

* (1 + eyi) (1 - £ 2 (yjyk + yjyt + ykye) - 2e 3 y j y k y e .) ' 

where (i,j,k,£) is any permutation of (1,2,3,4). 

It is not immediately clear that the map defined by formula (|4"U|) (or (J4"ll ) is birational. 
However, this is the case, and can be shown with the help of the following result. 

Theorem 8. For the map ^JD\j , there holds: 

R(y,e) ' 



R{y, -e) 
Proof. From (!4Up it follows that 



N N 



*=i i + fc »* i=i 1 fc y* 
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that is 



This is the claim of the Theorem. ■ 
Theorem [8] assures that discrete equations (jlOjl are equivalent to 

Vi = ^ - — ^ — -, 1 < i < N, 
1 + eyi 1 - eyi R(y, -e) 

which is nothing but O40p with y <-y y and e -H- — e. This means that the inverse map for 
f(y,s) is given by f~ 1 (y,e) = f(y, —e), and, in particular, that / is birational. 

The following statement shows that the map (14T1) admits the same conserved quantities 
and invariant volume form of the HK-discretization (|26|) . 

Theorem 9. 

(1) The functions 

Kij yi- yj y k yi 



K ke yiyj yk-yi 



l<i,j,k,£<N, (42) 



are conserved quantities of the map (f^7| ). Among these, there are N — 2 function- 
ally independent ones. 

(2) The map [Jlty admits an invariant volume form 

_ dyi A ■ ■ ■ A dy N 

where ip(y) is any one of the functions 

N-l / N ^ 2 



^) n 



J/* 

fe=i 



Proof. (1) It follows immediately from (j4*0j) that 



so that the quantity 



-~=r= 

y* % \3/i %/ 

y^-yy i/ii/j =J?( | (43) 
y<% y« - y? 

is independent of indices i,j. Therefore, functions (142]) are integrals of motion. 

(2) We have to compute the determinant det(dy/dy). To this end, we represent the 
map y i — y 'y as a composition of three simpler maps, 

y i — y u i — ^ v i — ^ y", 

where 



N - - 1 



Ui I 1 — £ > 
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One easily computes for the first and the third map: 

det (!) = rj(TrW' det (l)=rj(TTkp=n(i-* 

As for the remaining second map, we compute: 

N 



2 



Thus, 



3 3=1 



det (l^) = det [ ^(Rl + sue") ) = -J- ( l + c > 



duj \R 2K J R N \ R 7 R N + V 

Collecting everything and using (flO]) and (Ji"3"j) . we find: 

det UJ " Wia+^) 2 LU 2 



11 (1 + £yk) 2 11 ^2 

2/i - \ [rr^ 



n 



In the case N = 3, Theorem [2] tells us that map f l4T]) with admits the following conserved 
quantities: 

K mn (e) = K mn e), e) = , 

1 - ^ViVj 

where K mn are defined in ([3]) and i, j — 1, 2, 3, i ^ j. 
In the case iV = 4, we have the following result. 

Theorem 10. The map ( f77| ) w^A iV = 4 admits the following conserved quantities: 

K mn {e) = K mn e), $(y, e) = 1 (44) 

V(l - £ - £ 

ifere i^ mn are defined in ^EEj) , and (i,j,k,£) is any permutation of (1,2,3,4). Among 
functions [44\ )> there are three functionally independent ones. 

Proof. We have to prove that 

" y m ~ Vn VmVn \ 2 mMMh _ (1 - £ 2 M/j)(l - £ 2 ym) ^ 
y m y n Vm-Vn) 2/i2/22/32/4 (1 - 6 2 yiyj)(l - 6 2 y k y t ) ' 

For the l.h.s. of f l4"5]) . equations fj4Tj) and f T43|) immediately yield the following representa- 
tion: 

' Vm-Vn VraVn \ VtSzVzVA D 2 (y,6) 

(46) 



VmVn Vm-Vn J VmVzVi Di(y , s) D 2 (y , e)D 3 (y , e) D 4 (y , e)' 

Turning to the r.h.s. of fj4"5j) . we use Eq. (1411) for N = 4 to obtain: 



1 - 6 2 Viyj = 1 - 



e 2 ViVj{l + ey k ) 2 (l + eyt) 2 



D i (y,e)D j (y,e) 

Now, a direct computation proves the following polynomial identity (valid only if N 

4): 

Di(y, e)Dj(y, 6) - 6 2 y iyj (l + 6y k ) 2 (l + 6y e ) 2 = (1 - 6 2 y k y e )D(y, e). 
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Therefore, we find: 

1 - e^iVj 

1 - £ 2 ykVe 
so that _ _ _ _ 

(1 -e 2 y i y j )(l -e 2 y k ye) _ 

(1 -£ 2 ^yj)(l - e 2 y k y e ) Di(y,e)Dj(y,s)D k (y,s)D e (y,e)' 
Comparing (|46|) and (|47|) . we see that (|45|) holds true. ■ 

7. Concluding remarks 

There are several particular problems left open in the above exposition. The most 
important and difficult one is about the existence or non-existence of the last, (N — 
l)-st, integral of motion of our maps. However, still more important is to understand 
general reasons for integrability or non-integrability of birational maps generated within 
the Hirota-Kimura approach. A still more important probram is the study of integrability 
of birational maps on a possibly general basis. We hope that additional "experimental" 
results presented in this note will be helpful for this ambitious goal. 

Acknowledgements 

We would like to thank Vsevolod Adler for calling our attention to the generalized 
Kovalevskaya system and for communicating his observation that the cross-ratios of the 
phase variables are integrals of motion of the HK-discretization of this system. 

The authors are partly supported by DFG (Deutsche Forschungsgemeinschaft)in the 
frame of Sonderforschungsbereich/Transregio 109 "Discretization in Geometry and Dy- 
namics" . 

The participation of the second author in the SIDE 10 conference in Ningbo (China) 
was financially supported by the organisers, which is gratefully acknowledged. 

References 

[I] A.V. Borisov, I.S. Mamaev, Poisson structures and Lie-algebras in Hamiltonian mechanics, Izhevsk, 
1999. 

[2] D.B. Fairlie, An elegant integrable system, Phys. Lett. A, 119/9, 438-440, 1987. 
[3] R. Hirota, K. Kimura, Discretization of the Euler top, Jour. Phys. Soc. Japan, 69, 627-630, 2000. 
[4] A.N.W. Hone, M. Petrera, Three-dimensional discrete systems of Hirota-Kimura type and deformed 
Lie-Poisson algebras, Jour. Gcom. Mech., 1/1, 55-85, 2009. 
[5] Correspondence of S. V. Kovalevskaya and G. Mittag-Leffler, Nauka, 1984. 

[6] M. Petrera, A. Pfadler, Yu.B. Suris, On integrability of Hirota-Kimura type discretizations. Experi- 
mental study of the discrete Clebsch system, Exp. Math., 18/2, 223-247, 2009. 

[7] M. Petrera, A. Pfadler, Yu.B. Suris, On integrability of Hirota-Kimura type discretizations, Reg. 
Chaot. Dyn., 16/3-4, 245-289, 2011. 

[8] M. Petrera, Yu.B. Suris, On the Hamiltonian structure of Hirota-Kimura discretization of the Euler 
top, Math. Nachr., 283/11, 1654-1663, 2011. 

[9] M. Petrera, Yu.B. Suris, Spherical geometry and integrable systems, in preparation. 

[10] A.G. Reyman, M.A. Semenov-Tian-Shansky, Group theoretical methods in the theory of finite- 
dimensional integrable systems, in: Dynamical Systems VII, Springer, Berlin, 1994. 

[II] Yu.B. Suris, The problem of integrable discretization: Hamiltonian approach, Progress in Mathe- 
matics, Vol. 219., Basel: Birkhauser, 2003. 



D(y,e) 
D i (y,£)D j (y,ey 

D 2 (y,e) 



